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Abstract 

We compute the effects of soft-gluon resummation, at the next-to-leading-logarithmic 
level, in the fixed-target hadroproduction cross section for prompt photons. We find in 
general that the corrections to the fixed next-to-leading-order results are large for large 
renormalization scales, and small for small scales. This leads to a significant reduction of the 
scale dependence of the results for most experimental configurations of interest. We compare 
our results to the recent measurements by the E706 and UA6 collaborations. 
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1 Introduction 



The phenomenological interest of prompt-photon production in fixed-target experiments |T|, Q 
resides mainly in its use as a gluon probe in structure-function studies. Prompt-photon production 
is historically our main source of information on the gluon parton density at large x (e.g. x > 
0.2) [|[ H, ^ H, a region which has very little infiuence on the evolution of the deep-inelastic- 
scattering structure functions. This same region is relevant for hadron colliders in production 
phenomena at very large transverse momenta, and thus its understanding is crucial in order to 
disentangle possible new physics signals from the QCD background. 

For example, a particularly interesting problem has emerged in the past few years in the pro- 
duction of large-transverse-energy [Ex) jets at the Tevatron. An excess over the QCD prediction 
has been reported by the CDF collaboration 0, for jets with Ex > 350 GeV. While the excess has 
not been confirmed by the D0 data it is of interest to study the uncertainty in the high- £^7- 
tail of the jet distribution due to the gluon density systematics, to see whether there is room 
for deviations as large as those detected by CDF. For example, a suitable modification of the 
gluon density at large x has been proposed (CTEQ4HJ Jill), which is consistent with the excess 



observed by CDF. The study of the recent E706 prompt-photon data [g], however, suggests that 
a consistent fit of the large-x^ {xt = 2Et/VS) rate is incompatible with the CTEQ4HJ gluon 
density 0. Moreover, both jet cross sections and direct-photon cross sections at high transverse 
energy are affected by soft-gluon effects. These effects should be understood in both cases in order 
to be able to claim a discrepancy with QCD predictions. In particular, these effects can be very 
important in the direct-photon case, since the typical E-r values probed are much smaller than in 
the case of jet production at the Tevatron and, therefore, the size of the running coupling at 
the relevant scales is bigger. 

Comparisons between theory and prompt-photon experimental results have been carried out 
recently in Refs. |^, ^, 0, 13|. The recent E706 data 0] seem to differ most from the next- 



to-leading order calculation, over the whole xt range. In Refs. 0, [T^, an attempt is made 
to fit the E706 data by introducing an intrinsic transverse momentum of the incoming partons 
with ^1.2-1.4 GeV^ The precise details of how the intrinsic kx is incorporated in the 
calculations, however, can significantly affect the impact of these corrections, as shown by the 
large variations reported in Ref. [0. 

The use of an intrinsic-transverse-momentum model is sometimes motivated as a way of es- 
timating the effects of soft-gluon emission. The most prominent effect of soft-gluon emission in 
Drell-Yan pair production is the generation of the characteristic transverse-momentum spectrum 
of the lepton pair. This can be modeled with an appropriate intrinsic transverse momentum of the 
incoming partons. As a matter of fact, the formalism for soft-gluon resummation in Drell-Yan pair 
production can be shown to merge, at very small transverse momenta, into some non-perturbative 



intrinsic transverse momentum of the partons inside the hadron |T^. While this approach is not 
unreasonable when one considers the transverse momentum of the produced pair, it can however 
lead to inconsistencies for the problem of sm^f/e-photon production. In fact, for example, it is 
quite clear that the photon xt spectrum at large xt explores the kinematic region of x — 1 in 
the parton densities, which is certainly not the case for the transverse-momentum distribution of 
a Drell-Yan pair. Thus, as of now, a method for the inclusion of non-perturbative effects in the 
resummed formulae for the high-xy limit of the inclusive photon cross section is not available. 
Furthermore, in the opposite limit of small xt, it is the multiple emission of hard (rather than 
soft) gluons that leads to a sizeable perturbative broadening of the transverse momenta of the 
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incoming partons O, [16 



In this work, we consider the effect of soft-gluon resummation in prompt-photon production 
near the threshold hmit, that is to say for a;^ — >■ 1. The theoretical evaluation of these effects, at 
the next-to- leading logarithmic accuracy, has been carried out independently in Refs. fl^ and [|l^ . 
We shall review in the next section the necessary formalism, using the language of Ref. |TB|. In 
the rest of the paper, we will present its phenomenological applications, and we will thus discuss 
its numerical implementation as well as its impact on physical cross sections. 

As is well known, prompt-photon production takes place both by hard-photon emission from 
initial- or final-state quarks (direct component), and by collinear radiation from final-state partons. 
This last mechanism is not fully calculable in perturbation theory and, in fact, it depends upon the 
photon fragmentation function. Because of the large suppression of the fragmentation function at 
large momentum fractions 2;, it is usually believed that this contribution becomes irrelevant when 
xt increases. Contrary to common wisdom, we shall instead show that the very-large-x^ behaviour 
of the direct and of the fragmentation production processes is the same if the incoming hadrons 
do not contain valence antiquarks, as in the case of pN collisions. Under these circumstances, 
resummation should therefore be performed for the fragmentation emission too. We will show, 
however, that in the cases of practical interest the corrections due to the fragmentation processes 
are small, and we shall limit our considerations to the hard-photon part. 

The plan of the paper is as follows. In Section |^ we review the formalism for the resummation 
of threshold effects, and the main formulas valid for the specific case of prompt-photon produc- 
tion. There we also recall the main issues related to the inversion of the resummed expressions 
from Mellin space back to the physical x space. In Section |^ we study numerically the impact of 
the resummation corrections. We explore the effects both at the parton and hadron level, con- 
sidering kinematical configurations and distributions of phenomenological relevance for current 
experiments. In particular, we concentrate on the study of the size of the resummation correc- 
tions, and of the residual dependence on the choice of renormalization and factorization scales. 
Section |^ contains a comparison between our results and the data from some recent experiments. 
This does not want to be a comprehensive phenomenological study, but a preliminary analysis of 
the impact of our results on the comparison of theory and data. Our conclusions, and the outlook 
for future progress, are given in Section ^ An Appendix collects some details of the resummation 
formulas. 



2 Theoretical framework and notation 
2.1 Kinematics and cross section 

We consider the inclusive production of a single prompt photon in hadron collisions: 

Hi{Pi) + H2{P2) ^ lip) + X . (1) 

The colliding hadrons Hi and H2 carry momenta and P2, respectively. We parametrize the 
momenta in terms of light-cone coordinates: 

= (P+,Pt,P") , P^ = ^{P°±P^). (2) 
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In their centre-of-mass frame, using massless kinematics, the momenta of the colhding hadrons 
have the following light-cone coordinates 



Pi 



^ -1,0,0), P,^ = J^ {0,0,1] 



(3) 



where S = {Pi + ^2)^ is the centre-of-mass energy squared. The photon momentum p is thus 
parametrized as 

(4) 



P 



where and y are the transverse energy and the rapidity, respectively. We also introduce the 
customary scaling variable xt (0 < < 1): 

2Et 



Xt 



(5) 



In the present paper we are mostly interested in the prompt-photon production cross section 
integrated over y at fixed Ft- According to perturbative QCD, the cross section is given by the 
following factorization formula 



da^{xT, Et) 



z53 / dXi fa/Hr{Xl,liD / A/Ha (3^2, /^f) 

„u Jo Jo 



dx < 6 \ X 



Xt 



^XiX2 

+ y dz dc/^{z,fi'j) — 



Crab^-y{x, as{fi^);E^, /i^, fip, fij] 



(6) 



Xt 



Z^XiX2 



aab^c{x, E^, A^J) 



where a,b,c denotes the parton indices (a = q,q,g), and fa/Hi{xi, fip) and fb/H2{xi^ f^'p) the 
parton densities of the colliding hadrons, evaluated at the factorization scale fip- The first and 
the second term in the curly bracket on the right-hand side of Eq. represent the direct and 
the fragmentation component of the cross section, respectively. The fragmentation component 
involves the parton fragmentation function dc/'y{z, fij-) of the observed photon at the factorization 
scale fif, which, in general, differs from the scale fip of the parton densities. 

The rescalec^ partonic cross sections aah-.-y and aah^c in Eq- (§) are computable in QCD 
perturbation theory as power series expansions in the running coupling as{^i'^), /i being the renor- 
malization scale in the MS renormalization scheme: 



aab^^{x, as{fi ); Er^, fi , fip, fij) = a as{fi 



aab^c{x, a.(Ai^); E^, /ij) = 



n=l 



(7) 



n=l 



Note that the ratio between the direct and the fragmentation terms in Eqs. (0) and (|^) is of the 
order of a /as, where a is the fine structure constant. This ratio is compensated by the photon- 
fragmentation function dc/^y, which (at least formally) is of the order of a/ug, so that direct and 
fragmentation components equally contribute to Eq. (|^). 

■^These functions are related to the partonic differential cross sections by ^ab^i = Ej^ daab^i / dEr {i = 7,c). 
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Figure 1: Behaviour of the LO and NLO terms o'^^.^^^ and o-^^i^^^ (see Eq. (0)) of the 
direct component of the prompt-photon cross-section. The contributions of the partonic 
channels ah = qq (left) and ah = qg (right) are rescaled by the factor (1 — x^) and plotted 
as a function of x. The renormalization, factorization and fragmentation scales are all set 
equal to ji^ = 2E^, and = 1. 



Throughout the paper we always use parton densities and parton fragmentation functions as 
defined in the MS factorization scheme. In general, we consider different values for the renormal- 
ization and factorization scales /i, hf, l^f-, although we always assume that all of them are of the 
order of the photon transverse energy Ex- 

The LO terms cr^ab-^d-y -^1- ® following parton-scattering subprocesses at the 

tree-level 

q + q-* g + 1 ^ q + g ^ q + i , q + g-^q + i- (9) 

Using our normalization, the two independent (non- vanishing) partonic cross sections for the direct 
component are: 



- (0) / \ 2^F 



2 



(2 - x') (10) 



= = - ^ (l + , (11) 

where is the quark electric charge. Note that, having integrated over the photon pseudorapidity, 
the expressions (0) and ([TI|) are even functions of the photon transverse energy E^, i.e. they 
depend on rather than on x. The NLO terms o'^ab-^-y -'^l- (0) were first computed in Ref. [ pO[| . 

The partonic contributions ^ab^c to the fragmentation component of the cross section are 
exactly equal to those of the single-hadron inclusive distribution. Note that, unlike in the case 
of the direct component, all the parton-parton scattering subprocesses ah ^ c (i.e. including 
ah = qq,gg) contribute to the fragmentation component already at LO. The explicit calculation 



of CTab^c up to NLO was performed in Ref. |21 



The behaviour of the LO and NLO perturbative contributions to the direct component of the 
prompt-photon cross section is shown in Fig. 0. 
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The LO terms a'^^^^g^{x) and aq^J-^qj{x) are both singular when a; — >• 1: 



a^^_^,.^{x) -j=== , (x^l) , (12) 



and they both vanish in the high-energy hmit x ^ 0. The integrable singularity in Eq. ( p!2D is 
a typical phase-space effect, while the vanishing behaviour at small-x is due to the dominance of 
fermion (i.e. spin 1/2) exchange in the t-channel. 

Two new dynamical features appear at NLO. Near the threshold region x —>■ 1, the NLO 
contributions are double-logarithmically enhanced, 



cr 



^^\x) a'^''\x) \n\l - x) , (X-.1) , (13) 



because the radiation of soft and, possibly, collinear partons is strongly inhibited by the kinematics. 



In the high-energy limit x —>■ 0, the partonic cross sections o"^],!^^(x) approach constant values [ p!6[] : 
this Regge plateau follows from the fact that at NLO single-gluon (i.e. spin 1) t-channel exchange 
affects all the partonic subprocesses. The behaviour of the partonic contributions (3"^),^^(x) in the 

remaining intermediate region of x has no straightforward physical interpretation (e.g. aq]^-^-y{x) 
even becomes negative) because it strongly depends on the scale-dependent corrections already 
subtracted in the definition of the parton densities and parton fragmentation functions. 

Higher-order perturbative QCD corrections in the small-x^ regime can systematically be com- 



puted by using the fc^-factorization approach [15|, which consistently takes into account the per- 



turbative broadening of the transverse momenta of the incoming partons. 

We are interested in this work in the behaviour of the QCD corrections near the partonic- 
threshold region a; ^ 1, i.e. when the transverse energy of the photon approaches the partonic 
centre-of-mass energy y/xiX2S. In this region, the singularities in Eqs. (|1^, |13]) are enhanced 
by double-logarithmic corrections due to soft-gluon radiation and the higher-order cross section 
contributions in Eqs. ^ behave as 



a 



(")(x) ~ a(°)(a;) K.2n ln'"(l -x)+ a„,2„-iln2"-i(l - x) + . . .] . (14) 



Resummation of these soft-gluon effects to all orders in perturbation theory can be important to 
improve the reliability of the QCD predictions. 



2.2 A^- moment space 



The resummation program of soft-gluon contributions has to be carried out ^ in the 

Mellin-transform space, or A^-space. Working in A^-space, we can disentangle the soft-gluon ef- 
fects in the parton densities from those in the partonic cross section and we can straightforwardly 
implement and factorize the kinematic constraints of energy and longitudinal-momentum conser- 
vation. 



The latter point is particularly relevant for soft-gluon resummation in hadron collisions [0. 
Indeed, all-order soft-momentum recoil cannot exactly be taken into account by directly working 
in x-space and the ensuing kinematics approximation leads to (same-sign) factorially growing co- 



efficients. This implies [0 that no resummed logarithmic hierarchy can consistently be defined in 
x-space (the classes of leading logs ^^"(1 — x), next-to-leading logs ln^"~^(l — x) and so forth in 
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Eq. ( p!4D are not separately summable, because they lead to divergent and not integrable contri- 
butions at a; = 1). On the contrary, no kinematics approximation (in the soft limit) is required in 
A^-space and the corresponding logarithmic hierarchy of In iV-contributions is systematically well 
defined. 

To work in A^-space, it is convenient to consider the Mellin transform cr^ jv(-E'T) of the dimen- 
sionless hadronic distribution Ej, da.y{xT, Ex) / dEx- The A^-moments with respect to x"^ and at 
fixed Et are thus defined as follows: 



Jo 



dxi {xir-^ ^3^M^i^ _ (15) 

dEx 

In A^-moment space, Eq. (|^) takes a simple factorized form 

a,b 

X ja-afe^^, A^(a,(/i^);^|,/i^/i|,/iJ) (16) 

+ ^ (^ab^c, N{0is{fJ^'^)] E^, fl^ , fj,p, f/j) cic/7,2Af+3(/^/) 
c 

where we have introduced the customary A^-moments fa/H,N and da/^^^N of the parton densities 
and parton fragmentation functions: 

fa/H,N{f^^) = / dx X^~^ fa/H{x,fi^) , (17) 







da/-f,N{l^) = dz Z da/^{z,IJ,) . {U 







Note that the A^-moments of the partonic cross sections in Eq. ([T6|) are again defined with respect 



to x"^: 



o-ab^-y, Ar(a4/i^); E^, i/, /jj) = / dx"^ (x^)^ ^ aab^-f{x, asii^^); E^, /x^) 

^0 



(19) 



The explicit expressions of the A^-moments CF^qq^g^ ^qg->q'y n contributions in Eqs. (jT^) 

and (|lT]) were obtained in Ref. [|l^ and are recalled in Appendix 0. 



Note also the pattern of moment indices in the various factors of Eq. (|TBp, that is, fa/H,N+i 
for the parton densities and c?c/7, 2Ar+3 for the parton fragmentation functions. This non-trivial 
pattern follows from the conservation of the longitudinal and transverse momenta. 

The threshold region — *• 1 corresponds to the limit A^ — >■ cxd in A^-moment space. In this 
limit, the soft-gluon corrections (|14D to the higher-order contributions of the partonic cross sections 
become 

a!;^ ~ [ Cn,2n In'" N + Cn,2n-1 In'""^ A^ + . . .] . (20) 

The resummation of the soft-gluon logarithmic corrections to all orders in perturbation theory 
has been considered in Refs. and |T^. In the following section we recall the main results. 
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2.3 Soft-gluon resummation at high 



In Ref. [|T8| soft-gluon resummation has been performed in detail for the various partonic channels 
that contribute to the direct component of the prompt-photon cross section jv (Et) in Eq. 

We discuss first the large- iV behaviour of the partonic cross sections CTab^-y, n for the partonic 
channels ab = qq, qg, qg that start to contribute at LO. These cross sections can be written as 



^ (res) 



[1 + 0{as/N)] , ab = qq, qg, qg 



(21) 



where 0{as/N) denotes terms that contribute beyond LO and are furthermore suppressed by a rel- 
ative factor 0{1/N) at large N. The logarithmically-enhanced soft-gluon corrections are included 
in the resummed expressions d''^b^^, n t)e factorized with respect to the corresponding LO 

cross sections 0"^°^^^ ^. The a//- order resummation formulae are 



a 



(res) 



[as{n'^);El,fi^,n]„^)) 



a asijJ, ) a 



2\ -(0) ^ 



2\ 1,2. r\2 



a,{^^^),Q^/^^^■Q^/^^l) 



(22) 



a 



(res) 

qg^-y, N 



;a,(/i2);E|,/i^/i|,/xJ) 



a «s(/i^) ^fJ^q-y,N Cgg-.^/{as{n^),Qyfi^;Q^/fil) 



(23) 



a. 



(res) 



;a,(/i^);£:^,/i^/i^,/xj) = a 



(res) 

qg^"f, N 



;a,(/i^);£;^,/i^/i^,/ij) 



(24) 



where 



(25) 



The functions Cab-^^yicts) in Eqs. (|2^ , |23| ) do not depend on N. Thus, the In A^-dependence of 
the resummed cross sections is entirely embodied by the radiative factors A'^^'^'^ . They depend on 
the flavour of the QCD partons a, b, d involved in the LO hard-scattering subprocess a + 6 — > (i + 7 
and can be expressed in an exponential form: 



(q2 Q2\ ^ ^ 

"^(/i^), — ; — ) = exp i In N g'^^J {boa sifi'^) InN) 



+ 9abiboas{^J:') InN, Q'^/fi''; Q'' / ij^'p) + 0(«,(«, IniV) 



(26) 



where 60 is the first coefficient of the QCD /3-function 

IICa-^TrN 



f 



(27) 



Note that the functions g^^\ g^"^^ and so forth in the exponent do not depend separately on 
and IniV. They are functions of the expansion variable A = boaglTiN and vanish when A = 0. 
This means that the exponentiation structure in Eq. (^) is not trivial and, in particular, that 
all the double logarithmic (DL) terms a'^Cn,2n In^" in Eq. (^) are taken into account by simply 
exponentiating the lowest-order contribution a<(Ci_2 In^ N. The exponentiation in Eq. (26) defines 



an improved perturbative expansion in the threshold region. The function InN g^^^ resums all the 
/eac/mg' logarithmic (LL) contributions a"ln"'''^ in the exponent, g^"^^ contains the next-to-leading 
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logarithmic (NLL) terms a" In" A^, and so forth. Once the functions g^'^^ have been computed, 
we have a systematic perturbative treatment of the region of N where In ^ 1, which is much 
larger than the domain In^ <^ 1 where the fixed-order calculation in is reliable. 

The LL and NLL functions g^^'^ and g'^'^^ in Eq. (p6D have been explicitly computed in Ref. [|18 
The LL functions g^^'^ are different for the qq and qg partonic channels of Eqs. 



and ( P3D but 

they can be expressed in terms of parton colour factors and a single (parton-independent) function 



{2Cf - Ca) h^^\X) + Ca h^^\\/2) 
CAh^'\\) + CFh^'\\/2) , 



with 



/i«(A) 



1 



2nbnX 



2A + (1 -2A)ln(l -2A) 



(28) 



(29) 



The explicit expressions of the NLL functions g'^'^^ are recalled in Appendix 

Note that the LL functions g^^^ do not depend on the factorization scale fxp- This dependence 
starts to appear only in the NLL functions g^'^\ Note also the mismatch between the moment 
index of the radiative factor and that of o"l°L^„, m in Eqs. (|22|, p3l): the former depends on A^ + 1, 



ab^d'y, N 

like the parton densities in Eq. (0). The explicit /lir-dependence of g^'^^ exactly matches the 
scale dependence of the parton densities at large values of A^. Thus, when (and only when) NLL 



resummation is included, we can expect |18, 25| better stabilization of the calculation of the cross 
section at large xt with respect to variations of the factorization scale (see Sec. ^.2[ ). 

The functions Cab-^-yic^s) in Eqs. (p2|, |23|) contain all the terms that are constant in the large- A^ 



limit. They are produced by hard virtual contributions and by subdominant (non- logarithmic) 
soft corrections to the LO hard-scattering subprocesses. These functions are computable as power 
series expansions in 



n=l 



1+ ^C^^lJQV^^';QV^^l) + Oial) 



TT 



(30) 



gq^^ and CilL-y in Eqs. (H, H, H. 



At present, we know only the first-order constant coefficients C^^l 
These coefficients can be extracted |T^ from the complete NLO analytic results of Refs. pO| , p6| , p7[| . 
Their values are recalled in Appendix 0. 

The inclusion of the A^-independent function Cab^-yic^s) in the resummed formulae does not 
affect the shape of the cross section near threshold, but improves the soft-gluon resummation by 
fixing the overall (perturbative) normalization of the logarithmic radiative factor. 

We can explicitly show |]18|, ESl 



the theoretical improvement that is obtained by combin- 



ing the NLL radiative factor with the first-order coefficient C^),^^. Expanding the resummation 
formulae (|2^, |23|) in towers of logarithmic contributions as in Eq. (^Up , we have 



^ (res) / ,-, 



2 2 2 \ 



J 1 + 5^ «H C„,2„ In^- AT + Cn,2n-l{El/^ll, 



n=l 

+Cn,2n-2{El/ ^ll, / ^l^) lu^^'^ A^ + ^(In^"-^ | 



(31) 
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where = ^^(/i^). The dominant and next-to-dominant coefficients Cn^2n and 

Cn.2ra— 1 ^^^6 Con- 
trolled by evaluating the radiative factor to NLL accuracy. When the NLL radiative factor is 
supplemented with the coefficient C^|,L,^, we can correctly control also the coefficients c„,2n-2- In 



particular, we can predict [18| the large- behaviour of the next-to-next-to-leading order (NNLO) 



cross sections cr^^^^ in Eq. (|^) up to C(lnA^). 

Note also that the coefficients c„,2n are scale independent and the coefficients c„,2n-i depend 
on the sole factorization scale fip. In the tower expansion (|3TD , the first terms that explicitly 
depend on the renormalization scale fj, (and on fip, as well) are those controlled by Cn,2n-2- Their 
dependence on /x is obtained by combining that of C^^^^(2£'|>//i|., 2i?|>//i^) with that of the 

radiative factor at NLL order. The inclusion of the ffist-order constant coefficient C^),!^^ thus 
theoretically stabilizes the resummed partonic cross section at large xt with respect to variations 
of the renormalization scale. This scale dependence is numerically studied in Sec. p.2| . 

So far we have only considered the near-threshold behaviour of the partonic cross sections 



(Tgg^-y, AT, CTgg^^^N, o'qg^'^.N in Eq- (|2T|). The behaviour of other partonic channels ah ^ ^ that 



contribute to the direct component of the prompt-photon cross section was discussed in Ref. [p^ . 
It turns out that the partonic channel ah = gg enters the resummed cross section only at next- 
to-next-to-leading logarithmic (NNLL) accuracy and that all the other channels are relatively 
suppressed in the same way as the correction 0{as/N) on the right-hand side of Eq. (^if ). Since 
we are interested in explicitly performing soft-gluon resummation up to NLL order, we can limit 
ourselves to considering the resummed expressions in Eqs. (pTD -(p4D. 

Detailed numerical studies of the resummed cross sections are presented in Sec. |. However, 
from the analytical results reviewed in this section, we may already anticipate that soft-gluon 
resummation increases the perturbative QCD predictions in the large-xr region. This conclusion 
can be argued by a simplified treatment within the DL approximation. To DL accuracy, the 
exponent of the radiative factors in Eq. (EBI) has to be expanded to its ffist order in as, and we 



obtain 






^ (res) 
%^7, N 

-(0) 


~ exp < 




« (res) 
'^qq^-/, N 
^(0) 

^qq^qf,N 


~ exp < 




^ (res) 
^qg^-y, N 


^ 19-* 




^ (res) 
^gg^7, N 


- ^(0) 
qq^ 


<?7, 



\^ln'N]=exp{{CF + 2CA)^ln'N]>l , (32) 



2^ 



In^ at} = exp |(4Cf - Ca) ^ In^ A^} > 1 , (33) 



JO) 



exp {3{Ca - CF)^\n' n} > -fp^^^ . (34) 



"qq^qi, N 

For the sake of completeness, in the square bracket on the right-hand side of Eqs. (0) and (p3D we 
have explicitly separated the positive contributions coming from the initial-state partons and the 
negative contribution from the final-state recoil. From these equations we see that the resummed 
partonic cross sections o"^^'^-^, n ^"^^ ^qg^-y, n ^^^^ both enhanced with respect to their LO approxi- 
mations o'g^Lg-y AT) ^qg^q'y N' Moreover, the enhancement in the qg partonic channel is larger than 



that in the qq channel. We refer the reader to Ref. |]I8| for a discussion on the physical origin of 
this behaviour. 
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2.4 Fragmentation component 



We can now comment on the large- behaviour of the fragmentation component of the prompt- 
photon cross section, by comparing the direct and fragmentation contributions in Eq. ([161) . 

The partonic cross sections (jab^-y, n and aah^c, n have the same large- behaviour, but, owing 
to the hard (although collinear) emission always involved in any splitting process c — > 7 + X, the 
photon fragmentation function dc/y^N is of the order of 1/N. Therefore, in the curly bracket on 
the right-hand side of Eq. (|16[) the fragmentation component is formally suppressed by a factor 
of 1/N with respect to the direct component. This suppression is consistent with the fact that 



the resummed partonic cross sections for the direct processes (see the right-hand side of Eqs. (22) 
and (P^D) turn out to be independent of the photon fragmentation scale /i/. 

This argument shows that, in many cases, the fragmentation contributions are subdominant 
near threshold and, thus, they can be neglected in resummed calculations at large xt- 

The caveat 'in many cases' in the above conclusion regards the fact that the argument applies 
to the partonic contributions in the curly bracket of Eq. (|16D. In other words, the argument 
assumes that all the different initial-state partonic channels ab give comparable contributions to 
the hadronic cross section. This is not always true once the effect of the parton densities is 
included. 

A relevant exception is indeed the case of prompt-photon production in proton-nucleon colli- 
sions. Owing to the low antiquark content of the colliding hadrons, the hadronic cross section is 
mostly due to the partonic channels ab = qg and ab = qq: 

a^,^(ET) ~a^;^^(^T) + <,v(^T) , (piV collisions) . (35) 

As for the qg initial-state contribution a'^^j^{ET), we can use the above argument to conclude that 
its direct component dominates at large Et- Setting all the scales equal to Ft, for the sake of 
simplicity, we can write: 

<;v(^t) - <rH^T) - U,N-,i{El) Um^i{E'^) N{as{E'^)) ■ (36) 

However, in the case of the qq initial state, the direct component enters only at NLO and, thus, the 
cross section is dominated by the fragmentation part and, in particular, by photon fragmentation 
from a final-state quark of the LO scattering subprocess q + q ^ q + q. We can write: 

a'^'^iEr) ~ <r'^(i^;T) ~ f,,N+iiE^) kN+i{El) a,,^,, nMED) d,,,^2N+,{El) . (37) 

Taking the ratio of the two initial-state contributions and replacing &qq~>q^ n{(^s) and ^qg^-y, Af (o^s) 
by their LO contributions in Eqs. (^, H), we obtain 

The factor ag/cx on the right-hand side is compensated by the behaviour of the photon fragmen- 
tation function dqj^ 2Af+3 oc a/ as- In the large- limit, the ratio of the LO partonic contributions 
is constant and, thus, the fragmentation function produces an (9(l/A^)-suppression factor. 
Nonetheless, this suppression can be balanced by the parton density contribution fq,N+i/ fg,N+i 
since, at large x, the gluon density is typically softer than the quark density. As a matter of 
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fact, using the Altarelli-Parisi evolution equation at LO and under reasonable assumptions on the 
large-x behaviour of the parton densities at the initial evolution scale, it is easy to show that we 
have the following asymptotic behaviour at very large values of and of the evolution scale: 

ivliV ^^'^ ' ^^^^ 
la , ^ 

N inJy as 

Combining these results with Eq. (p8D, in the large- iV limit we therefore get 

, — constant . (41) 

This discussion shows that, in the case of large-xr prompt-photon production in pN colli- 
sions, the contribution of the fragmentation component to the hadronic cross section can become 
comparable to that of the direct component. 

Although on the right-hand side of Eq. (^) we have approximated the partonic contribu- 
tion d'qq^g^ ]\f{as) /d-qg^^^ ]\f{as) by its LO expansion, the inclusion of higher-order terms and, in 
particular, of resummation effects does not substantially modify the conclusion. 

Performing soft-gluon resummation in the partonic cross sections CTab-^c, n of the fragmentation 
component, we can write an expression that is analogous to Eq. (|2T|): 

^ab^c,N = ^i'XA^ + 0{ajN)] . (42) 

Limiting our treatment to the LL accuracy, the resummed cross section is given by [ll8| , |17 

The radiative factor is 

^-^-'^4as{^i'),%■%X] =exp{lnAr£)_^,^(6oa.(/i^)lnAr) + 0((a,lnAr)'=)| , (44) 

where the LL function g^^b^dc analogous to those in Eq. (^) and it can be expressed in terms 
of the colour charges Ca {Ca = Cp, if the parton a is a quark, and Ca = Ca, if a is a gluon) of the 
partons involved in the LO hard-scattering subprocess: 



9. 



(1) 

ab^dc 



( A) = (a + a + - Q) {X)+C,h^'\X/2) . (45) 



In particular, for the q + q ^ q + q channel we have 

9iU,W = h^'^ (A) + Cp h^'^ (A/2) , (46) 

which is very similar to gi^J{X) in Eq. ( p8| ) because 2Cf — Ca- More precisely, since 2Cp = 

Ca{1 — l/^c); 9(i^q-^ggW is slightly smaller than gq^JiX) as long as they are evaluated in the 
perturbative region A = boas^TiN < 1/2. 
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We can now come back to the effect of the fragmentation component in pN colhsions. Using 
the resummed partonic cross sections in Eqs. (p3| , |43|) rather than their LO approximations, the 
right-hand side of Eq. (|38| ) has to be multiphed by an additional contribution, as given by the 
ratio of the corresponding radiative factors, namely 

A^Vr ^ exp{(C^ + 4C^)tln^iV} 

A?.Tr"exp{(C^ + 2C^)tln^iV} " ^ 

Because of the relation ACp ^ between the colour charges, this factor does not sizeably 
differ from unity, as it can be argued by its DL approximation on the right-hand side of Eq. (|47| ) 
(see also the comment below Eq. (^61)). We have thus shown that, at least at the LL level, soft- 
gluon resummation does not enhance the relative importance of the fragmentation component for 
prompt-photon production at large Xt- 

The importance of the fragmentation component in pN collisions mainly depends on the de- 
tailed behaviour of the parton densities at large x and on how large are the values of of interest. 
This issue, as well as the impact of the NLL corrections to the LL results obtained above, require 
further studies that will be presented in a future work. As for the present study, we limit ourselves 
to perform soft-gluon resummation in the direct component and we check that the fragmentation 
component does not sizeably contribute to the hadronic cross section in the actual experimental 
configurations investigated in the paper (see Sec. 



2.5 Resummed cross section to NLL accuracy 

We use soft-gluon resummation to NLL accuracy at the parton level to introduce an improved 
prompt-photon cross section a^^^{ET) as follows 



a. 



(res) 
7,iV 



(Et) - ^ fa/Hi,N+l{f^F) fb/H2,N+l{f^F) 

ab=qq, qg ,qg 



a 



(res) 



>,(/i2);E|,/i^/x|) - (^ai;!5^^^(a,(/i2);E|,/i^/i|,)j^ 



+ a. 



(NLO) 

7, AT 



{Ej 



(48) 



where cr^jf^ is the prompt-photon hadronic cross section at NLO, cr^bll^ tv given in Eqs. (p^)- 
(|2^) and [a^^^fl^ ^ j represents its perturbative truncation at order aa'^ (i.e. at NLO). Thus, 



because of the subtraction in the square bracket on the right-hand side, Eq. (|48| ) exactly reproduces 
the NLO results and resums soft-gluon effects beyond O(aa^) to NLL accuracy. In general, we 
evaluate o"^^^^^ n using the NLL expression ( p6D of the radiative factors and including the 0{as) 
contribution ( PU] ) of the constant factors Cab^-y{cis)- This defines our NLO-I-NLL predictions. 

The resummed formulae presented so far are given in A^-moment space. To obtain cross sections 
in the physical xy-space (i.e. as functions of the centre-of-mass energy), one has to perform the 
inverse Mellin transformation: 



3 da!y''"\xT,ET) 
^ dEr 



1 
2Tii 



Cmp+«oo 



dN x; 



-2Af^{res) 



{E_ 



T 



(49) 
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When the A^-moments are evaluated at a fixed perturbative order in they are analytic 
functions in a right half-plane of the complex variable N . In this case, the constant Cmp that 
defines the integration contour in Eq. (^9|) has to be chosen in this half-plane, that is, on the right 
of all the possible singularities of the A^-moments. 

An additional complication occurs when the A^-moments are computed in resummed pertur- 
bation theory. In this case, since the resummed functions g^ab^^) -^'l- O) (^^ NLL 
functions g^^^) are singular at A = 1/2, the soft-gluon factors AAr(Q;s(/i^)) in Eq. (^) have cut 
singularities that start at the branch-point A^ = A^^^ = exp(l/2feoas)- These singularities, which 
are related to the divergent behaviour of the perturbative running coupling as near the Landau 
pole, signal the onset of non-perturbative phenomena at very large values of A^ or, equivalently, 
in the region very close to threshold. 

The issue of how to deal with the Landau singularity in soft-gluon resummation formulae for 
hadronic collisions was discussed in detail in Ref. [jl9[. In the evaluation of the inverse Mellin 



transformation (^) we thus use the Minimal Prescription introduced in Ref. [|r^. The constant 
Cmp is chosen in such a way that all singularities in the integrand are to the left of the inte- 
gration contour, except for the Landau singularity at A^ = Ni, that should lie to the far right. 
This prescription is consistent |T^ with the perturbative content of the soft-gluon resummation 
formulae because it converges asymptotically to the perturbative series and it does not introduce 
(unjustified) power corrections of non-perturbative origin. These corrections are certainly present 
in physical cross-sections, but their effect is not expected to be sizeable as long as Et is sufficiently 
perturbative and xt is sufficiently far from the hadronic threshold. Obviously, approaching the 
essentially non-perturbative regime ~ 1 GeV, physically motivated treatment 

of non-perturbative effects has to be introduced. In the following sections, we limit ourselves 
to presenting numerical and phenomenological results that do not include any non-perturbative 
correction. 



3 Results 

We present in this section some numerical results, to provide an illustration of the size of the 
effects considered and to show the improvements obtained with respect to scale variations after 
the inclusion of the NLL corrections. 



3.1 Parton-level results 

We start by discussing the resummation effects at the level of partonic cross sections. The re- 
summed partonic cross section can be obtained from Eqs. (^) and ( ^9]) by assuming parton-density 
functions of the form fa/nix) = 5(1 — x), and hence fa/H,N = 1 for all complex values of A^. 

We consider first the 0{aa1) terms in the expansion of the resummed cross section, in order to 
estimate to which accuracy this reproduces the exact NLO results. In Fig. || (left) we plot the func- 
tion a^qq^^/E'^, defined in Eq. (^, as a function of 77 = (1 —xt)/xt- The exact 0{aa1) result pO| 
is compared with three possible implementations of the resummation procedure, all equivalent at 
NLL accuracy. The first case (short-dashed line) corresponds to our default resummed prediction, 
as given by Eq. (^21). In the second case (dot-dashed line), we set the constant C^q-,^, introduced 
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Figure 2: Left (Right): the 0{aa1) contribution to the partonic prompt-photon cross 
section for the process qq j + X {qg — > 'j + X) plotted as a function of r/ = {1 — xt)/xt- 
The sohd hne represents the exact NLO result of Ref. pO]; the short-dashed line is the 
0{aa1) piece of of the resummed result defined by Eqs. (|22|) and (p3|); the dot-dashed line 
is obtained from this last result by setting the constant Cg^l^^ (Cqll^'y) to 0; the dashed 
line is obtained using Eq. (^) (Eq. (|5ll)), with A = 2. The renormalization, factorization 
and fragmentation scales were all set equal to fi'^ = 2E^, and Cq = 1. 



in Eq. (|30D equal to 0. In the third case, we keep the contribution of the constant Cg^l^^, but 
we modify it by a term suppressed by a factor of in order to explore the possible effect of 
contributions of order which cannot be taken into account by the soft-gluon resummation. 
As a constraint on the form of these corrections, we must impose that no poles appear on the 
positive real axis in the N plane (these poles would logarithmically enhance the partonic cross 
section when xt — 0). We select a parametrization of the corrections that allows us to 
bracket the exact result at (9(aa^) : 

^^S-l. - CS^, (1 + ^^^_, ) , (50) 



- Cill^ (1 - , (51) 

with y4 > 0. In our applications we shall consider the two cases with A = (namely no correction 
to the Cab^'y term) and A = 2 as a way to establish the size of subleading threshold corrections 
beyond the NLL order. 

As one can see from Fig. ^, the inclusion of the finite term Cg^^^ is essential to accurately 
reproduce shape and normalization of the exact 0{aa^g) result not only near threshold, but below 
it as well. The agreement deteriorates unavoidably for 77 ^ 1, as, here, terms subleading in 
become important. 

Analogous results for the qg channel are given in the right panel of Fig. ^ Note that in both 
cases the two choices A = and A = 2 m Eq. (|5T|) bracket the exact result over a large region 
of ?7, and thus provide a good estimator of the subleading terms' systematics. The choice A = 0, 
furthermore, provides a very accurate description up to values of of the order of 1/10. 



14 



400 



200 



-200 
10" 



- 


1 


"1 1 1 1 ' 
&qq (pb/GeV) 




1 Et=10 


GeV //-\ 
//- \ 
/■/ \ 


- 




/J \ 

'A 

" / "\ 






/ dots: LO ^ 
/ dashes: NLO 
^ solid: NLO+NLL (A=0) 

dotdash: NLO+NLL (A=2) \ 




1 





'"■' Y I I I I I 

- dotdash: NLO+NLL (A=2) (l-x^) 5-,g (pb/GeV) 

100 — solid: NLO + NLL (A=0) 

: dashes: NLO 

ao — dots: LO 




?7 = (1-Xt)/xt 



?] = (1-Xt)/xi. 



Figure 3: Partonic cross-section for the processes qq "fX (left) and qg — > 7X (riglit) 
(in pb/GeV, and for = 10 GeV). The dotted hne is the LO result; the dashed line 
is the exact NLO result; the solid (dotdashed) lines correspond to the NLO+NLL result, 



with the coefficient A defined in Eq. ( |50D (left) and in Eq. (51) (right) equal to (2). The 
number of flavours was set equal to 4 and we have taken Aq^^d = 0.151 GeV. 




7] = (1-Xt)/Xt 77 = (1-Xt)/xt 

Figure 4: Same as Fig. |, for Et = 100 GeV. 
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Et (GeV) Et (GeV) 

Figure 5: Contribution of gluon resummation at order 0{aa^) and higher, relative 
to the exact NLO result, for photon production via qq (left plot) and qg (right plot) 
annihilation, in pN collisions at ^/S = 31.5 GeV (-Ebeam = 530 GeV). The solid (dashed) 
lines correspond to ^4 = (A = 2). The three sets of curves correspond to the choice of 
scale fi = hf = '^Et, Et and Et/2, in descending order, with PDF set CTEQ4M and 
Nf = 5. 

The fully-resummed parton-level cross sections are shown in Fig. ^ for the qq and qg channels 
(left and right panel, respectively), and for Et = 10 GeV (Fig. ^ collects the same results for 
Et = 100 GeV). Here and in the following we shall define the resummed cross sections as in 
Eq. (^5D, that is, we substitute their Ol^aa"^) terms with the exact NLO result, using the same 
choice of renormalization and factorization scales. In this way our results are exact up to (and 
including) 0{aal), and include the NLL resummation of terms of 0{aal) and higher. We compare 
the fixed-order results (dashed lines) with the resummed results. For these we provide both the 
A = and A = 2 prescriptions. Note that, even at the level of resummed cross sections, the 
difference between the A = and A = 2 results are rather small, in particular for the qq channel. 



3.2 Hadron-level results 

In this section we present some results for the full hadronic cross sections. The main points we 
intend to highlight are: 

1. the size of the NLL corrections, relative to the NLO contributions; 

2. the scale dependence at NLL order. 

Our goal here is to explore the pure effects of resummation at higher orders. Therefore we 
shall neglect in this section all production channels which are not improved by the resummation 
corrections considered in this work. This includes all processes which first appear at 0{aa1), 
such as gg — > qq'j and qq^'^ as well as all contributions proportional to a parton 7 

fragmentation function, as discussed in the previous section. These terms will however be included, 
at fixed NLO, in our comparison with experimental data, performed in the next section. 

To be more specific, we list here the classes of diagrams included at NLO, in addition to the 
LO processes qq —>■ (77 and qg — > q'j (the possible replacement of quarks with antiquarks in all 
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Figure 6: Same as Fig. |5|, for the combined production channels qg + qq. 



these cases is understood): 



qg (191 , m mi , m 



Hi 



m 991 



(52) 



For the third set we only include the diagrams proportional to e^, since the part of the amplitude 
describing the photon emission from the final-state quarks (which is by itself gauge invariant), 
cannot be considered as a correction to any tree-level process. 



As a default set of parton densities we shall use the CTEQ4M set described in Ref. For 
the purposes of the present study, no significant change is obtained if different sets are used. 

Figures ^ and ^ present the ratios: 



NLL, qq 



(T. 



NLO 
qq 



0" 



NLL, qg 



0", 



NLO 



^res _ ^NLO 

'^NLL, (qg+qq) ^(qg+qq) 



NLO 

qq 



cr, 



NLO 

qg 



cr, 



NLO 

(qg+qq) 



(53) 



where, for simplicity, we indicated here with a the differential distribution dajdEj^. For each 
channel we present the results using both the A = (solid lines) and A = 2 (dashed lines) 
prescriptions. We also show the dependence on the choice of renormalization and factorization 
scales, which we take equal, and varying within the set fi = fip = {Et/2, Et, ^Et). In this section 
we shall always keep the fragmentation scale Hf, necessary for the factorization of the singularities 
from final-state coUinear photon emission, equal to Et- Note that the size of the resummation 
effects is larger for the larger scales, contrary to the behaviour of the scale dependence of the 
NLO cross section. This suggests that the scale dependence of the resummed cross section will be 
reduced relative to that of the NLO results. 

The scale dependence of the resummed cross section [qg + qq contributions), compared to the 
NLO one, is given in Fig. |^, for the A = case. The same result for A = 2 is given in the left panel 
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Figure 7: Scale dependence of da/dEx {qg + qq components) for prompt photons in 
pN collisions, at Et = 5 GeV, plotted as a function of the proton-beam energy, -Ebeam 
(the associated values of xt are given on the top scale). The solid lines represent the 
exact NLO result for different choices of = fj-p (fJ- = Et/2 and 2Et), normalized to 
the fi = Et result. The dashed lines represent the NLO+NLL result (with A = 0) for 
different choices of /x, normalized to the NLO fj. = Et result. 
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Figure 8: Same as Fig. |^ but with A = 2 (left panel) and with cj^il^^ 



(right panel). 
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Table 1: Contributions to the prompt-photon rate da/dEx in pN colhsions at -Ebcam = 
530 GeV, from higher orders in the expansion of the NLL resummed resuh. Resuhs 
for Et = 5 (10) GeV are shown in the first (second) two rows. The renormahzation, 
factorization and fragmentation scales are set equal to /U = = /^/ = Et, and the PDF 
set is CTEQ4M. The aa^ column gives the exact NLO result. 



of Fig. ||. We plot the distributions as a function of the beam energy (-Eboam) for the fixed value 
of Et = 5 GeV. Different values of -Ebeam, therefore, probe different ranges of xt, as indicated by 
the upper labels on the plots. 

Note the significant reduction in scale dependence, more marked in the A = case. To display 
the importance of the inclusion of the constant C^],!^^ terms, we show the same scale-dependence 

plot with C^fel^^ = in the right panel of Fig. ^ While the scale sensitivity is slightly worse than 

in the cases with C^J ^ 0, there is still an important improvement over the NLO behaviour. 



Similar results, for Et = 10 GeV, are shown in Figs. and 10. The general features of these 



distributions are similar to those of the plots for Et = 5 GeV. Small violations of xy-scaling can 
be observed between Et = 5 and 10 GeV, due to the evolution of the coupling constant and of 
parton densities. 

We also explored the independent renormahzation- and factorization-scale dependence of our 
calculations. The large size of this dependence at NLO was stressed already in Refs. 13|. The 



results, for pN coHisions at -Ebeam = 5 and 10 GeV, are shown in Figs. |TT] and jT^], respectively. 
With the exception of the renormalization-scale dependence at 5 GeV, a significant improvement 
in the stability of the results, relative to the dependence at LO and NLO, is observed in all cases. 

The convergence of the higher-order corrections is displayed in Table |I|. The last column 
includes the sum of all contributions of order aa^ and higher, performed using the Minimal 
Prescription of Eq. (^9]). The fixed-order terms do not have any ambiguity due to the choice 
of the contour for the Mellin transformation in Eq. (^9]) . The contribution from the qq channel 
converges very rapidly. In the case of the qg channel the convergence is slower, in particular at 
the larger values of xt, but even at Et = 10 GeV the size of the resummed contributions beyond 
order aa^ is only of the order of 10% of the total. This supports the validity of the Minimal 
Prescription, since the truncated resummed expansion converges to it very smoothly. 
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Figure 9: Scale dependence of da/dEx {qg + qq components) for prompt photons in 
pN collisions, at Et = 10 GeV, plotted as a function of the proton-beam energy, E'beam 
(the associated values of xt are given on the top scale). The solid lines represent the 
exact NLO result for different choices oi (j, = jip (/i = E^j^ and 2Et), normalized to 
the fi = Et result. The dashed lines represent the NLO+NLL result (with A = 0) for 
different choices of normalized to the NLO /i = Et result. 
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Figure 10: Same as Fig. but with with A = 2 (left panel) and with C^^^^ 
panel). 
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Figure 11: Scale dependence of the differential Et distribution in pN collisions, for 
Et = 5 GeV. We compare the results at the Born and NLO level with the results of 
the resummed calculation. Upper left: renormalization-scale dependence, with the fac- 
torization scale fixed to /Ltp = Et- Upper right: factorization-scale dependence, with the 
rcnormalization scale fixed to ^ = Ex- Lower left: scale dependence, with the renormal- 
ization and factorization scales equal. 
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4 Comparison with current data 



We present in this section a comparison between our calculations and the results of the two most 
recent measurements of fixed-target direct-photon production. E706 studied photons produced 
in pBe collisions at -Ebeam = 530 GeV, covering the centre-of-mass rapidity range \y\ < 0.75f\ 
The data span the Et region between 3.5 and 11 GeV, approximately. This corresponds to 
0.22 < xt < 0.70. UA6 studied photons produced in pp collisions at -Ebeam = 315 GeV, over the 
rapidity range —0.1 < y < 0.9. The Et region extends between 4 and 7 GeV (0.33 < xt < 0.58), 
approximately. What follows is not meant to be a systematic phenomenological study, but will 
serve as a benchmark to assess the impact of the resummation effects in realistic experimental 
conditions. All the calculations in this section have been done assuming /x = /ip = /i/. For 
recent complete studies of all available data, done using the fixed-order NL calculations, see 
Refs. 0, §. 

To compare our resummed predictions with actual data, two additional things need to be done: 
inclusion of the 1/A^-suppressed contributions, and inclusion of realistic experimental cuts. 

The class of processes for which we evaluated resummation corrections in section ^ provides 
the dominant contribution to the production rate in realistic experimental configurations. On the 
left-hand side of Fig. |T3| we plot the relative contribution, evaluated at NLO, of the processes with 
a qg, qq, qq^'"* and gg initial state for the E706 experimental configuration. As one can see, the sum 
of qg and qq accounts for 90% of the overall rate, independently of Et- On the right-hand side of 
Fig. 13 we plot the rate of the direct contributions, relative to the sum of direct and fragmentation, 
for each given channel. In the case of the gg^'-* and gg channels we compare the absolute values 
of the rates, since the direct component is negative after the subtraction of the initial-state mass 
singularities. The comparison of the two plots in Fig. |l^ shows that the processes for which we are 
going to include resummation corrections account at NLO for a fraction of the total rate between 



70 and 90%, in the Et range 4-12 GeV. The situation is even better in pp collisions (see Fig. ^ 



obtained for the UA6 experimental configuration), where the to-be-resummed processes account 
for over 90% of the NLO rate. 

We therefore expect that the neglect of the resummation corrections to the gg + qq^'^ and to the 
fragmentation processes is only a minor correction to the overall picture. For the present study, and 
in addition to the resummed predictions for the processes listed in Eq. (^2]), we will therefore only 
include the fixed-order NL determination of these remaining components of the direct-photon 
production process. The fragmentation processes are evaluated using the NLO single-inclusive 



parton Et distributions from |2T| , convoluted with the GRV photon fragmentation functions 
We found very small sensitivity to the choice of the photon fragmentation functions. 

As anticipated above, the comparison of our results with actual data requires the inclusion 
of realistic detector acceptance cuts. The resummation formalism discussed so far allows the 
evaluation of the transverse-energy distributions integrated over the full range of rapidity for 
the observed photon. This approximation is technically correct, provided the measurement is 
performed within a finite range in rapidity including the value of ?/ = in the collision centre- 
of-mass frame. This is because in the large-xr limit all production is concentrated at y = 0. To 
include the effect of experimental rapidity cuts, which usually do include the y = point, we 



^In this section we shall use y to indicate the value of photon rapidity in the hadron-hadron centre-of-mass 
frame. 
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Figure 13: Left: relative contribution of different initial states to the Et distribution 
in pN collisions at E'bcam = 530 GeV. Right: relative size of the direct contribution vs. 
the sum of direct and fragmentation one, for the different channels. 
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Figure 14: Left: relative contribution of different initial states to the distribution in 
pp collisions at £?beam = 315 GeV. Right: relative size of the direct contribution vs. the 
sum of direct and fragmentation one, for the different channels. 
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Figure 15: Effect of the rapidity cuts, as a function of Et- Upper set of curves: E706. 
Lower set of curves: UA6. The three curves in each set are obtained with [i = E^j^^ Et 
and 2Et. 



therefore apply the following acceptance correction to our resummed cross sections: 

(NLO),' Y] 

(yeY)^ (all ,) ^(,,ol|ln ,) " (^4) 
The experimental configuration for most experiments of practical interest is wide enough that, for 



the relevant values of Et, the rapidity acceptance is very large. We show two examples in Fig. ^ 

One set of curves gives the ratio <j{\y\ < 0.75)/cr(all y), evaluated at NLO for the E706 
experimental configuration. The three curves correspond to different choices of scale n, and show 
very small dependence on fi. We also checked that the dependence on the PDF set used is at the 
level of 1-2%. The acceptance loss is of the order of 25% for Et values around 4 GeV, and becomes 
totally negligible for Et > 8 GeV. The other set gives the ratio cr(— 0.1 < y < 0.9)/cr(all y), 
evaluated at NLO for the UA6 experimental configuration. The acceptance loss is here more 
significant, due to the tight cut at negative rapidity. 



We present our prediction for the E706 data in Fig. |16p|. Notice the significant reduction in 



scale dependence obtained when going from the fixed-order NL calculation to the resummed result. 
This scale reduction is particularly evident at high Et, where the resummation effects are more 
important. Notice that while at low Et the band with the resummed prediction is all contained 
within the NLO uncertainty band, at high Et the NLL result becomes larger relative to NLO for 
all the displayed scale choices. The plot shows a reasonable agreement between data and theory 
at large Et, indicating that no additional significant contribution is required in this region. The 
large disagreement between data and theory already present at NLO is still present, as no net 
increase is obtained from the resummation contributions. Their only effect is to reduce the scale 
dependence. 



^The theoretical prediction for pN has been rescaled by a 1.09 factor, to account for nuclear corrections to the 
pBe process M. 
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Figure 16: E706 data compared to the resummed theoretical predictions. The theory 
was rescaled by a factor of 1.09 to account for nuclear corrections in the conversion from 
the pBe to the pN rate For comparison, the figure includes as well the fixed-order 
NL result. We used PDF set CTEQ4M, and GRV photon fragmentation functions. 



A similar picture emerges from the comparison of the theory with the pp UA6 data. This is 



shown in Fig. 0. The disagreement between data and theory at low Et is however much less 
dramatic here than in the case of E706. The extent to which these Iow-Et discrepancies can 
be removed by the inclusion of non-perturbative effects such as an intrinsic remains to be 
understood, as the global consistency of the different data sets is not very compelling O, O B[. 



5 Discussion and conclusions 



We presented in this paper a numerical study of the impact of resummation corrections, at the 
next-to-leading logarithmic level, on the transverse energy distribution of direct photons produced 
in hadronic collisions. We dealt with the resummation of the — >■ 1 Sudakov logarithms studied 
theoretically in Refs. |18, 17]. As a result, this work is mostly of relevance for typical fixed-target 
photon production. The current prompt-photon data from the high-energy hadronic colliders 
cover in fact the region xt ^ 0.1, where Sudakov effects are negligible. 

We showed that the inclusion of higher-order Sudakov corrections improves significantly the 
factorization and renormalization scale dependence, relative to what observed in the fixed-order 
NL calculations. Even when the scales are varied independently, the uncertainty from scale vari- 
ations is significantly reduced. As a result of the reduced scale dependence, the overall size of 
the resummation contributions depends significantly on the chosen scale. In general, however, 
the resummed cross sections for different scale choices have values contained within the NLO 
uncertainty band for xt values up to 0.5, and exceed the upper side of the NLO band by large 
factors when xt approaches 1. Still, our resummation corrections turn out to be much smaller, at 
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Figure 17: UA6 pp data compared to the resummed theoretical predictions. For com- 
parison, the figure includes as well the fixed-order NL result. We used PDF set CTEQ4M, 
and GRV photon fragmentation functions. 

least in the range of existing data, than the effects induced at large xt by some implementations 
of intrinsic kx effects, as discussed in Refs. |0, |T^, 0. In these papers the effect of intrinsic-fcy 
corrections was evaluated to be as large as factors of 3 and more, for the whole Et range. We 
believe that this is related to the absence in the intrinsic-fcr models of the appropriate Sudakov 
suppression due to the presence of the hadronic system recoiling against the photon (represented 
at LL in our formalism by the negative terms in the exponents of Eqs. (|32D and (^)). 



In our work we did not include resummation corrections to the fragmentation processes. We 



qq 



(') 



qq(')j 



proved that, in pN collisions, the large- behaviour of the corrections to the 
processes is formally similar to that of the corrections to the leading one, qg — ^7. These 
corrections are therefore not suppressed when N increases towards larger values. Whether they 
can be neglected or not, is therefore a pure matter of numerics. We showed that their contribution 
is not dominant in the Et regions of experimental interest, and limited ourselves to including them 
at the fixed next-to-leading order. As discussed in Sec. |2.4|, we have no reasons to believe that the 



resummation corrections are any larger for the fragmentation processes than for the qg channel. 
A more quantitative study of these statements, and a complete phenomenological assessment of 
the comparison between theory and the current sets of data in view of the results presented in 
this paper, will be the subject of future work. 



A Appendix: Formulae for the resummed cross section 



In this Appendix we recall (see Ref. [^]) the explicit expressions of the various factors that 
contribute to the resummed cross sections in Eqs. ( |22D and ([23|). We use the customary notation 
for the colour factors in SU{N^) QCD, namely, Cp = (N^ - l)/(2Ar,), Ca = and Tr = 1/2. 
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The A^-moments of the LO partonic cross sections in Eqs. (|T0|, |TT]) are 



.(0) _ r(i/2) r(iv + i) 

^qq^g-y, N - ^ '=^q 



. (0) _ . (0) 

^qg^q-y^N ^ qg^q-y, N 



Nc T{N + 5/2) 

2 1 r(i/2)r(iv + i 



2 + N) , 

(7 + 5N) 



(55) 
(56) 



"SA^c r(A^ + 5/2) 
where T{z) is the Euler F-function. 

The formulae for the NLL functions g^'^^ in the exponent of the radiative factors in Eq. ( [26| ) 
are the following 



(2)/ 0^,0^ 



+ 



i2CF - Ca) Z^^'HA) + 2Ca h^^\\/2) 

^^^ln21n(l - 2A) + ^^^^ln(l - A) - ^ Aln^ 



(57) 



2A + ln(l - 2A) 



+ 



2TTh 



>- 



21n(l - A) - ln(l - 2A) 



In^ 



(:7^/i(2)(A)+2(:7^/i(2)(A/2) 



(58) 



+ — -In21n(l-2A) + — ln(l - A) Aln — 

ItTOo AtCOq TTOo jJ^p 



+ 



27cbo 



2A + ln(l-2A) 



+ 



Cp_ 
27r6o 



21n(l - A) - ln(l - 2A) 



In 



where 'jp = 0.5772 ... is the Euler number and bo, hi are the first two coefficients of the QCD 
/3-function 



, _ IICa - ATrN} _ VJC\ - IQCaTrNj - GCpTrNj 

The auxiliary function h^"^^ that appears in Eqs. (p7| , has the following expression 



(59) 



hi 
2Tihl 



2A + ln(l - 2A) + - ln2(l - 2A) 



where the coefficient K is given by 



K = C 



A 



67 



6 



10. 



2A + ln(l-2A) , (60) 



(61) 



The first-order coefficients CgJ^^ and Cqg-^-y of the A^-independent functions in Eq. (pOf ) are 



Cill,{QVfi';QVfil) = jl{2Cp-^CA)+lE[7rho-{2Cp-CA)\n2]-^i2Cp-CA)ln2 

+ ^K-Kg + ^(2Cp-^CA)+li2Cp-CA)\n''2 (62) 



- (2-fpCp - ^Cp] In % - nho In 
V 2 / fij, 
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-CF-CA\n2 



1 2 1 

- -K, + —(2CF + 19CA)+-CFln^2 

2 ^ 60\ ^) 2 



— (C^ -2Ca) In 2 
10^ ' 



- (7e{Cf + Ca) - ^Cf - vr6o) In % - nbo In ^ 



where 



2 6 ' 



(63) 



(64) 



and the coefficient K is given in Eq. (|6l]). 

Note that the LL functions g^^^ are given in Eq. (^). Thus, the formulae presented in this 
Appendix complete all the ingredients that are necessary to evaluate the resummed cross sections 
with NLL accuracy. 
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